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1. Introduction 

In this paper we study the third order nonselfadjoint spectral problem 

— <j> (x) — zm(x)(j)(x) for x G R, 

<t> x (-oc) = (f)xx(-oo) = 0, <t>xx{oo) = 0, 

which we refer to as the cubic string with Neumann-like boundary conditions. 
In particular we consider the case when m(x) is a discrete measure (a finite 
linear combination of Dirac deltas), and the problem of reconstructing m(x) given 
appropriate spectral data. The purpose of this introduction is to explain where this 
problem comes from and why it might be of interest. 

Our starting point is the Degasperis-Procesi (DP) equation, which is the 
completely integrable nonlinear PDE 

u t ~ u txx + iuu x = 3u x u xx + uu xxx . (1-2) 

It was discovered purely mathematically in a search for integrable equations [lj [2] , but 
has been shown by Johnson [3] to play a role in shallow water wave theory similar to 
that of the Camassa-Holm (CH) equation 0] 

u t — utxx + 3uu x = 2u x u xx + uu xxx . (1.3) 

By taking appropriate scaling limits one can derive new nonlinear equations which 
inherit some of the integrable structure of their ancestors. Consider in particular the 
high-frequency limit obtained by changing variables x i— > ex, t i— > st, and then letting 
e — > 0. For the CH equation (|1.3|) this gives the Hunter-Saxton equation [5J H>] for 
nematic liquid crystals, 

(u t + uu x ) xx = u x u xx . (1.4) 

The DP equation Q1.2p . on the other hand, reduces in the high-frequency limit to 

(u t + uu x ) xx = 0, (1.5) 

which we will refer to as the derivative Burgers equation. The name of course 
refers to the well known inviscid Burgers equation (or Riemann shock wave equation) 
u t + uu x = 0, the prototype equation for studying shock waves. That (| 1 . 5[) . with the 
two extra derivatives, has a Lax pair is a fairly recent discovery; as far as we are aware, 
this was first shown by Qiao and Li [7] (who studied the whole integrable hierarchy 
that (|1.5p belongs to). 

A unified way of writing these four equations is 

m t + m x u + bmu x = 0, (1-6) 

with m and b as in the following table: 



6 = 2 



b = 3 



Camassa-Holm 



Hunter-Saxton 



derivative Burgers 



Degasperis-Procesi 

(According to Dullin, Gottwald and Holm 8 , the whole family (|1.6[) with m — u — u xx 
and 6^—1 can be derived as an asymptotic limit of the shallow water equations. The 
CH and DP equations are the only integrable cases in this family. Moreover, (|1.6p with 
m = u — u xx is included in the more general family Ut + lSu+S'yu 2 /2+u 1 ~ u (v?u x ) x ] x + 
vu txx = studied by Rosenau [5]; take 1 + 2uj = b, —3j = 6+1, 8 — 0, v = — 1, 
and let x t— > —x. Another generalization of (II. 6|) is the pulson equation introduced by 
Fringer and Holm [TOl fTT] .) 
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A similar table summarizes the particular kinds of soliton-like solutions 
characteristic of these equations, as well as the spectral problems (including the type of 
boundary conditions) employed when computing these solutions using inverse spectral 
methods: 





Peakon solutions 


PlECEWISE LINEAR SOLUTIONS 


Discrete string 


Camassa-Holm 


Hunter-Saxton 


Discrete cubic string 


Degasperis-Procesi 


derivative Burgers 




Dirichlet(-like) 


Neumann( -like ) 



This requires perhaps some explanation. Peakons (peaked solitons) are solutions 
taking the form 



u(x,t) = 2jmfc(t) exp(- \x - x k (t)\), (1.7) 

k=l 

while the piecewise linear solutions are given by 

n 

u(x,t) = } m k (t) \x - x k (t)\ , (1.8) 



k=l 



where in each case the functions Xk(t) and rrik(t) (positions and momenta of the 
solitons) are required to satisfy a certain system of 2n ODEs in order for u(x, t) 
to satisfy the PDE. (See Theorem 12.31 below for the case of the derivative Burgers 
equation.) 

It was discovered in [12] that the spatial equation in the Lax pair of the CH 
equation can be transformed by a change of variables to the string equation with 
Dirichlet boundary conditions, 

— (j) (x) — zni(x)<f){x) for —1 < x < 1, 

#-i) = #i) = o, (L9) 

where z is the spectral parameter (squared frequency), m is the mass density and <f> 
represents the amplitude of the harmonic mode corresponding to to that frequency. 
(This originates of course in separation of variables in the linear wave equation 
v X x = rnv t t describing small vibrations of a string with given mass density m(x).) 
The peakon solutions are related to the case when the mass distribution is a sum 
of point masses (Dirac deltas); this is what is meant by the discrete string. Krein 
P31 HU [15] studied the inverse spectral problem of the string equation (| 1 .9|) and 
showed in particular that in the discrete case the mass distribution can be recovered 
from spectral data using continued fractions of Stieltjes type [16j . This was exploited 
in [T71 US] to solve the CH peakon ODEs explicitly and analyze the solutions in detail. 
Similarly, the inverse problem for the string equation on R with Neumann boundary 
conditions <p x (—oo) = 4> x (oo) = was used in [19] to solve the ODEs for the piecewise 
linear solutions of the Hunter-Saxton equation. The inverse spectral theory for the 
simple looking ODE (jl.9p is surprisingly rich in mathematical content. The most 
comprehensive account of the inverse problem for this equation, and also of Krcin's 
work on the subject, is the monograph [20] by Dym and McKean. 

Despite the superficial similarity, the DP equation has a mathematical structure 
quite different from that of the CH equation. In particular, the spatial equation in 
the Lax pair of the DP eqation is of third order instead of second. In [21] we showed 
how to transform it into the cubic string with Dirichlet- like boundary conditions, 

—6"'(x) = zm(x)6(x) for — 1 < x < 1, , 
y \ > (1.10) 

0(-l) = ^(-1) = 0, 0(1) = 0. 
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This problem was studied in great detail in our previous paper |22j , where the inverse 
problem was solved in the discrete case, leading to explicit solution formulas for the 
DP peakons. 

In this paper we first briefly show how the derivative Burgers equation is related 
to the cubic string with Neumann-like boundary conditions Then we will give 

a fairly complete account of the direct and inverse spectral theory of (jl.lll in the 
case when m(x) is a positive discrete measure. Some of the ideas used in [22) are 
developed further, and we also provide new shorter proofs for some key steps which 
perhaps appeared mysterious there. For the most part this paper is self-contained, 
with the exception of using some Heine-like multiple integral formulas proved in [22] 
(see Appendix B) and postponing the proof of Theorem 14.201 to another paper. 

2. The derivative Burgers equation 

Here we give a short account of the basic facts about the derivative Burgers equation 
(jl.5p . First we show how it inherits the integrable structure of the DP equation (jl.2p . 
The following Lax pair for the DP equation was given in [2] : 

(d x - d x )<p = zm<p, <j> t =[z~ 1 {c-dl)+u x -ud x \(t>, (2.1) 

where c is an arbitrary constant (or function of z); indeed, the compatibility condition 
'Ptxxx = <\>xxxt is equivalent to m t + m x u + 3mu x = and m x = (u — u xx ) x . The 
high-frequency substitution x h- > ex, 1 1— > et, together with m i— > — e~ 2 m, z i— > — e~ 1 z, 
and c i— > — ze~ 2 c, yields 

(e 2 d x - d x )4> = zrruf), <j> t = [z~ 1 dl +c + u x - ud x ] <j>, 

which is compatible iff mt + m x u + 3mu x — and m x = {u xx — e 2 u) x . Letting s — > 
we obtain the system 

- dl<j) = zm<p, (2.2) 

4> t = [z^d 2 x +c + u x - ud x ] 0, (2.3) 

which is compatible iff m t + m x ii + 3mu x = and m x — u xxx , and hence in particular 
when the derivative Burgers equation holds. 

Remark 2.1. Here the cubic string (|2.2p on the whole real line R appears naturally. 
This is in contrast to the DP case [21] where a change of variables mapping the real 
line to the finite interval (—1,1) is required to get rid of the d x term in the first 
equation of (|2.ip . A similar distinction holds for the ordinary string equation in the 
Hunter-Saxton and CH cases. 

Remark 2.2. Hone and Wang [23] have described a different high-frequency limititing 
procedure taking the DP equation to (the x derivative of) the integrable Vakhnenko 
equation (u t + uu x ) x + u — [24, 25J. It should perhaps also be mentioned that in the 
low-frequency limit e — > oo the substitution above directly reduces the DP equation 
to Ut + uu x = [2], but this does not give a useful Lax pair. 

We will return to the question of boundary conditions for (I2.2p shortly. The 
piecewise linear solutions of the derivative Burgers equation are described by the 
following theorem, where we use the convention sgnO = 0. Dots denote 4i as usual. 
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Theorem 2.3. The function u given by jl.8\) satisfies the derivative Burgers eguation 
\1.5\) in the sense of distributions if and only if 

n n 

Xk = y^ra» \%k — %i\ i mk — 2^ m k m l sgn(x 4 - x k ) (2.4) 

1=1 1=1 

for k — 1, . . . , n. 

One can assume that all m k ^ 0, since it follows from (|2.4p that any vanishing 
mfe remains zero and hence never enters the solution. 

The solutions in Theorem 12.31 satisfy the usual Burgers equation u t + uu x = 
(without the extra derivatives) only if X)fc=i m k = 0- This case will not be considered 
here, since we will make the following basic assumptions throughout this paper: 

X\ < ■ ■ ■ < x n , nik > for all k. (2.5) 

It can be shown that if these two assumptions hold for t — 0, then they will continue to 
hold, and the solution of (|2.4| will be defined, for all t > 0. (But if mt's of both signs 
are present, shocks may form after finite time, and then the solution will no longer 
have the form p.8[) .) Under the assumptions (|2.5p , sgn(iEi — Xk) can be replaced by 
sgn(i - k) in (j2T4"|) . 

Proposition 2.4. M = 53fe=i m k an d M+ = X)fe=i m kXk o,re constants of motion of 
the ODEs Pll ). 

Proof. This is easily verified by a short calculation. □ 

Since we are now in the context of piecewise linear solutions where u has the form 
(|1.8p . it follows that m = u xx = 2J^™ ra^ S(x — x^) is a (time-dependent) discrete 
measure. Consider some fixed time t. The wave function (f>(x, t; z) satisfying the cubic 
string equation (|2.2[) is piecewise a quadratic polynomial in x, since (j> xxx — away 
from the support of m. At the points Xk, (|2.2[) implies that </> and <fi x are continuous 
while <f) xx has a jump of size —2zmk <fi(xk)- These requirements completely define </> 
once </>, <p x , and <p xx are prescribed at some point. 

Now let u evolve according to the derivative Burgers equation; that is, let Xk and 
nik evolve according to (|2.4p . Then </> evolves according to (|2.3p . We can postulate 
that 

<t>(x,t;z) = 1 for x<x x {t) (2.6) 

provided that we choose c in (|2.3p equal to the constant of motion M = ^mfc; then 
(|2.6p is compatible with the time evolution since it makes both sides of (|2.3p vanish 
for a; < xi. Propagating (|2.6[) to the right gives 

4>(x, t; z) = A{t- z) ^ ~ + z) (x - x n ) + C(t; z) (2.7) 

for a; > x n (t), where A, B, C are some explicitly computable polynomials in z of 
degree n, with coefficients depending on the Xk(t)'s and mk(t)'s; see (13.31) below. 
Inserting this into (|2.3p . together with tt = mk(x — Xk) — Mx — M + for x > x n {t), 
gives the time evolution of these quantities: 

A 

A = 0, B = MB, C = — + IMC. (2.8) 

We see that for x > x n we can impose <p xx — 2 A = 0; this condition will be preserved 
by the evolution in time. Together with (|2.6p . this amounts to <f) satisfying the 
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Neumann-like boundary conditions <fi x (— oo) = (j> xx (— oo) = 0, <p xx (oo) = 0, which 
is how we were naturally led to the cubic string boundary value problem (jl.lj) . 

More details about the solutions of the ODEs (|2.4p will be given in a separate 
paper devoted to the derivative Burgers equation. Here we will leave this subject after 
showing that M is the first in a naturally appearing sequence of constants of motion, 
while M + is something of an odd bird. Indeed, writing A(z) = 2 ^^ =1 (— z) fc Mfc it is 
straightforward to derive an expression for the coefficients M k using (|3.3| . and since 
A is constant in time, we then obtain the following result. 

Theorem 2.5. The functions Mi, . . . , M n given by 

fc-i 

m *= e (ii m 0(nK -^ + j 2 )- (2.9) 

'e( [1 r') ieI 

are constants of motion for the ODEs ft2.4\ )- Here (^L™') denotes the set of all It- 
element subsets I = {ii < • • • < of the integer interval [l,n] — {1, . . . ,n}. The 
empty product when k = 1 is to be interpreted as 1, so that Mi = ^ TO fe = -W- 
For example, when n = 3 the constants of motion are 
Mi = + m 2 + m 3 = M, 

M 2 = mim 2 (xi - x 2 ) 2 + mim 3 (i:i - x 3 ) 2 + m 2 m 3 (x 2 - x 3 ) 2 , (2.10) 
M 3 = mim 2 m 3 (xi - x 2 ) 2 {x 2 - x 3 ) 2 , 
together with M+ = m\Xi + m 2 x 2 + TO3X3. 



3. The forward spectral problem 

We now turn to the spectral theory of the cubic string with Neumann-like boundary 
conditions The problem makes sense for any m(x) decaying sufficiently fast 

as I a; I — > 00, but we will deal exclusively with the discrete case where m(x) — 
2J2 m k $( x ~ x k), and in addition we assume ordering and positivity as in (|2.5[) . 
Then the eigenvalues z are the zeros of the time-independent polynomial A(z) defined 
by (|2.7j) . so the ODEs (|2.4[) associated to the derivative Burgers equation induce an 
isospectral deformation of If we suppress the t dependence in the notation, the 

propagation of the wave function from left to right is described explicitly by 

--G k {z)*(xi), (3.1) 



where $ 




1 
G k (z) =| 1 I . (3.2) 
-2m k z 1 





<$>{x+-z) = G n (z)L n _iG n -i{z)...L 2 G 2 {z)LiGi(z)\ . (3.3) 



From this one can easily extract the coefficients of A(z) = 2 Y^k=i(~ z ) k ^k explicitly; 
see Theorem [23] above. Note also that, as z — > 0, 
B(z) = 2z(M + - Mx n ) + 0{z 2 ), 



C(z) — 1 - z'^2 / m k {x n - x k ) 2 + 0(z 2 



(3.4) 



fc=i 
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Theorem 3.1. // m k > for all k, then the spectral problem hl.l}) with m(x) = 
2 ^2 m k 5(x — x k ) has n distinct nonnegative eigenvalues z — \ k , 

= A < Ai < ••• < A„_i. (3.5) 

Proof. It is trivial that z — is an eigenvalue, with corresponding eigenfunction 
0=1. The proof of the other assertions amounts to showing that the remaining 
eigenvalues coincide with the eigenvalues of a certain oscillatory matrix. The basic 
facts about oscillatory matrices are summarized in Appendix A. By (|3.ip . we have 
= G k+ i(z)L k ^(x+). With {C kl B k ,A k f = *(a£), this says that 

I 2 

Cfe+i = Cfe + l k B k + — -A k , 

B k +i = B k + l k A kl 

A k+1 = -2zm k+1 C k+1 + A k . 

The boundary conditions at — oo can be expressed as Aq = Bq = 0. The first equation 
gives 

2zC k+1 -2zC k = 2zl k B k + zl 2 k A k , 
which, using the second and third equations, becomes 

A k+ i - A k A k - A k -i ^-4 . 2 

1 = 2zl k } LA.+zLAk. 

m k+1 m k f-f 

The boundary condition at +oo is A — A n — 0, so X — {A\, , . . , ^4 n _i)' satisfies 
M.X — zCX, where M. is the symmetric tridiagonal matrix defined by A4u — 
m~ 1 + otLTj and Aij^-i = A^i-i^ = — m,^ 1 , while C is the lower triangular matrix 
with La = i? and dj = 2ljj for i > j. (A4 and C are both of size (n — 1) x (n — 1).) 
An easy induction on n shows that 

1 , . v mi + m 2 H h m„ 

detA4 = 



mim 2 . . .m„ 

which is positive when all > 0, so Mr 1 exists. Hence the nonzero eigenvalues of 
(jl.ip are the reciprocals of the eigenvalues of the matrix M~ 1 £. 

By Theorem IA.2I we are done if we can show that Ai' 1 C is oscillatory. The 
star operation does not alter the determinant, since Ai* = TAiT with T = 
diag(+l, —1, +1, . . . , ±1). Hence the determinant of Ai* is positive, and so are also 
all its principal minors, since they are determinants of the same form. Thus Ai* 
is positive definite, and Theorem IA.5I then implies that Ai* is oscillatory. Hence 
M' 1 = ((X*)*)" 1 is oscillatory by Theorem [Ol 

The matrix C is the weighted path matrix of the planar network in Figure [TJ 
and consequently totally nonnegative by Theorem IA.8I (Lindstrom's Lemma). Since 
C is also obviously nonsingular, Ai' 1 £ is oscillatory by Theorem IA.31 The proof is 
finished. □ 

Next, we define the Weyl functions W{z) and Z(z) of the discrete cubic string 
with Neumann-like boundary conditions, as well as the residues b k and c k in their 
partial fraction decompositions. The case M = ^2 m k = which would need special 
treatment is excluded here by our standing assumption that all m k > 0. Recall that 
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Figure 1. A planar network of order n — 1 whose weighted path matrix is C. 
Unlabelled edges have weight 1. 



the zeros of A(z) are z = and z = Afe, and that the lowest order terms of S(^) and 
C(z) are given by (|3.4|) . 

Definition 3.2. Let 



W(z) 



Z(*) 



^(a;n+) -6(2) 



E 



<f>(x n +) _ C{z) _ 1/2M 



(3.6) 



E^ 



fc=l 



z - Afe ' 



Theorem 3.3. The Weyl functions satisfy the relation 

Z(z) + Z{-z) = W[z)W(-z). (3.7) 

Hence, the second Weyl function Z is determined by the first Weyl function W ( except 
for the residue at z — which is determined by the constant M ) through the formula 

bjb k 



c k = W(-\ k )b k = -J2 



^ Xj + Afe 



k = 1, . . . , n — 1. 



Proof. Let 



S(z) = G n {z)L n - X G n -x{z) . . . L 2 G 2 (z)L 1 G 1 (z) 



(3.8) 



(3.9) 



denote the transition matrix from (|3.3|) . whose first column is (C, B,Ay. The (1,3) 
entry of the identity 



S(-z)*JS(z)J = I, 



where J 



1 
0-10 

1 



(3.10) 



reads 



A(-z)C(z) - B(-z)B(z) + C(-z)A(z) = 0. 



(Equation (|3.10p is proved by noting that X = G k (z) and X = L k both satisfy 
X(-z)- 1 = (JX(z)J)\ therefore so does X = S{z) 22.) By definition, W = B/A 
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and Z = C/A, so division by A{z)A(-z) yields (|3~7)) . The residue of (|3"77|) at z = A fc 
is (El- □ 



Remark 3.4. This proof (with trivial modifications) also works for proving the 
corresponding relation between the two Weyl functions for the the discrete cubic string 
with Dirichlet-like boundary conditions, which is used in the DP case. It is more 
straightforward than the proof in [22] , which relied on the analogue of Theorem 13.71 
below. 

Theorem 3.5. If rrik > for all k, then the residues bk of the Weyl function W(z) 
are all negative. 

Proof. By definition A, B, and C are all positive when z < 0. Hence W(— Afc) > 0, so 
that bk — B(Xk)/A'(\k) and Ck = C(Xk)/A'(Xk) are of the same sign by (|3. 8[) : they 
cannot be zero, since that would imply that A = B = C = for z — Afc, contradicting 
detS(z) = 1. 

We now use induction on n. For n — 2 the Weyl function is 

W( ) - B ^ ~ -2mili z -l/m 2 h h 

~ A(z) ~ -2(mi + m 2 )z + 2(m 1 m 2 ll)z 2 ~ _ mi + m 2 ~ z - X x ' 

Z mimsl? 

so b\ < 0. For n — p > 1, if we let mi — » + then W = W p tends to the corresponding 
Weyl function W p -± for the cubic string with parameters m 2l . . . , m p and l 2 , . . . , l p -\. 
For the largest eigenvalue and corresponding residue, we have A p _! — * +oo and 
6 p _! — > 0, while the other eigenvalues and residues tend to the eigenvalues and 
(negative, by induction hypothesis) residues of W p —\. Since the A^'s and b^s depend 
continuously on the parameters, the bk's cannot change sign during this process, and 
thus they must have been negative to begin with, except perhaps for b p -\. 

To prove that b p ^i < 0, note that B, by (|3 .3(1 and the definition of G p (z) and 
L p -i, does not depend on m p . By making m p small, we can make A p _i as large as 
desired; in particular, larger in absolute value than all zeros of B, which except for 
z = are also the zeros of W. 

Since the highest coefficients of A and B differ in sign, W < for all sufficiently 
large z > 0. If were positive, then W would be positive for some z slightly 
larger than A p _i, and that would force W to have a zero which is larger than A p _i, 
in contradiction to what we found above. Thus b p ^i < 0. □ 

Remark 3.6. Rather than using one spectrum of eigenvalues together with the 
corresponding residues, it is clear from (|3.6| that the Weyl function W(z) could instead 
be defined by two spectra, namely the zeros of 4>x{x n +) and <fixx(%n+)- Theorem 13.51 
above implies that these two sets of zeros are interlacing. In this regard, and also in 
being determined by a single Weyl function (which is perhaps surprising), the cubic 
string is similar to the classical Sturm-Liouville problem and to the ordinary string. 

Although it is not needed for the purposes of this paper, we also note how the 
Weyl functions change when the configuration {xk, mk}k=i of the discrete cubic string 
evolves in time according to the ODEs (|2.4p . 

Theorem 3.7. The time evolution of the Weyl functions is given by 

W(z) = MW(z), Z(z) = -+2MZ(z). (3.11) 

z 
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Hence, the residues evolve according to 

b k (t) = b k (0) e Mt , c k (t) = c fc (0) e 2Mt . (3.12) 

Proof. This follows immediately from (|2.8p . □ 



4. The inverse spectral problem 

A given discrete cubic string configuration {x kl m k } k l =1 with X\ < • • • < x n and all 
mj. > defines, as described in the previous section, a set of spectral data consisting 
of the positive eigenvalues Ai < • • • < A„_i together with their corresponding residues 
?)„_! < and the constant of motion M — J2 m k > 0. Only the relative 
positions l k = x k+ \ — x k enter this computation, which is why 2n — 1 is the correct 
number of spectral data if we wish to reconstruct the 2n — 1 string parameters 
h, . . . , l n -i and mi, ... , m n . 

Let a( m - ) denote the product of the first m factors in the transition matrix S{z) 
of (13~gi) : 

a {2k) (z) = G n (z)L n -iG n -i{z)L n - 2 ■ ■ ■ G„_ fe+ i(z)L„_ fc , 

(4.1) 

a {2k+l \z) = G n {z)L n _ x G n ^x{z)L n _ 2 . . .G n _ k+1 {z)L n _ k G n _ k (z). 

The Weyl functions W = B /A and Z = C/A are constructed from the first 
column <&(a:+) = (C, B,A) 1 in S(z) = a^ 2n ' 1 \ From the entries in the odd partial 
products a( 2fe+1 ) [z), which are polynomials in z, we will show how to construct 
rational approximations to W and Z with a common denominator, and satisfying 
certain normalization and symmetry properties rendering the approximants unique. 
Conversely, given a set of spectral data, we define W and Z by (|3.6|) . From the set of 
uniquely determined rational approximants referred to above we reconstruct all the 
factors Gj and Lj and thus solve the inverse problem. 

First, we have a simple power counting lemma. Recall that $ denotes the vector 
(4>, 4> x ,<j> xx ). For < k < n - 1, let 

Pk(z) \ 

Qk(z) =*(x-_ k ;z), (4.2) 
r k{z) ) 

so that (C,B,AY = §(x+) = a^ k+1 \p k ,q k ,r k y for any k. 

Lemma 4.1. If we compute the polynomial degree of a matrix entry-wise, then 

I k k- 1 k- 1 \ 
dega (2fc+1) (z) = k k-l k-l (4.3) 
\ k+l k k ) 

with the proviso that for k = the degrees in the second and third columns are all 
zero (negative degrees counted as zero). Furthermore, Aeg{p k ) = deg{q k ) = deg(r k ) = 
n — k — 1. 

Proof. Since each Gj(z) is linear in z, and since the leftmost factor G n (z) in (14. ip 
only affects the third row of a^ 2k \z), the first two rows of a^ 2k \z) have degree k — 1 
while the third row has degree k. Subsequent multiplication by G n - k (z) on the right 
gives the first claim. The second claim follows from an easy induction on k! = n — k 
which starts with (p n -i, Qn-i, r n -i) = (1,0,0) for k 1 = 1. □ 
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Fix k and write a(z) instead of cS 2k+1 \z) for simplicity. We now turn our 
attention to approximation properties of the entries of a(z). These properties will 
later model the approximation problems required to solve the inverse problem. 

Theorem 4.2. For k fixed, and z — > oo, 



a-31 



Z \ _ ( a n \ ( 0{z- v ) 

W ) \ 02! J V (!) 



Z \ ( Oi3 \ _ ( 0(z 



wher 



a 21 (0) ; ^ ; ' ^ a 22 (0) ; ^ 1 y ' (45) 

asi (0) = a 32 (0) = 0, a 33 (0) = 1. 
Moreover, setting W*(z) = — W{— z) and Z*(z) = Z{—£) we have 

a 3j Z* + a 2j W* + aij = 0(z-( fe+1 >), 1 < j < 3. (4.6) 

Proof. Regarding (|4.4p we only give the proof for the third column of a{z) to illustrate 
the main steps. Write (p,q,r) for (f>fc, (J^Tfe). Since 



we have 





031^23 - a 2 ia 33 p a 22 a 33 - a 32 a 23 q 



W _ ^21 — a 2lP + «2 2 g + a 2 3^ _ ^23 _ 033 r CL33 T 

a 3 3 a3iP + 032? + a 33 r a 33 a 3i 2 + a 32 f + a 33 

The 2x2 minors of a(z) appearing here are the (2, 1) and (1, 1) entries of the adjoint 
matrix adj(a(z)J, which equals a(z)~ 1 since each factor in a(z) has determinant 1. 
Hence the matrix 

adj(a(z)) = G n -k(-z)L~^ k G„-k+i(-z) . . . L~^ 2 G n -i(-z)L~l 1 G n (-z) 



has the same degree structure as in Lemma|4TTJ so the minors above both have degree k, 
and so do a 32 and a 33 . Moreover, p/r and q/r are both of order 0(1) by Lemma [4~T1 
while a 3 i has degree k + 1, so the whole expression is of order 0(z - ( fc+1 )) as claimed. 
The proof for Z — ai 3 /a 33 is similar, except that the appearing minors come from the 
first two rows of the second column of adj(a(z)) instead, which lowers the degree one 
step. 

Equation (I4.5[) follows from (|4.1[) since Gj(Q) — I. As for the symmetry relation 

(z-( fc+1 )) 
0(z~( fe +i)) 
(z-( fe +i)) 
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Like S(z), the matrix a(z) satisfies the property (|3.10p (the proof is the same), so the 
left-hand side above equals Ja(— z)* J (Z, W, 1)*, the rows of which give rise to (I4.6[) 
after changing z to — z. □ 

Proposition 4.3. For k = 0, 

a£\z) = -2m n z, a$(z)=Q, a$(z) = 1. (4.7) 
For k = 1, . . . , n — 1, 



"31 



71 

i—n—k / \i—n-\-l — k 



{2k+1) (z) = -2z( J2 H +--- + (-2z) k+1 ( J] ) "*»-*> (4.8) 



a 



(2fc+l) 
32 



(z) = -2z ^ mi(xi - x„_ fc ) 



Vi— n+1 — k 



+ • • • + (-2z)M Yl m i J ^- m n+1 _ fc Z„_ fe , (4.9) 



-i— n+2 — 

2 ^ 



4 2 3 fe+1 )(z) = l-2zf V m,. ( ^- X "- fc) 



»(*) = l-2*( E ... 2 

Vi=n+1 — k 



/ n i 2 \ 

+ ... + (-2z)M J] . (4.10) 

\i=n+l-fc / 

Proof. Write Gj = I — 2zm j (0 ) 0, 1)*(1, 0, 0) in the definition gj]) of a( 2fc+1 \ and 
collect the lowest and highest powers of z. □ 

Corollary 4.4. IFii/i [■] denoting the coefficient of the highest occurring power of z, 
m„_ fe = - 1 0<fc<n-l, (4.11) 

*"-» = g +1) ■ 1<*<»-1. (4-12) 

L«32 (*)J 

Next, we show that given W and the properties in Theorem 14.21 determine 
a(z) uniquely. Like in the DP case [22j we are dealing with a problem of simultaneous 
rational approximation, with approximants having a common denominator and for 
which the functions to be approximated are related by the quadratic constraint (|3.7p . 
In addition to the actual approximation conditions (14.41) . there is also the symmetry 
condition (|4.6p which makes the approximants unique despite the low order of the 
approximation condition. This type of approximation seems to be quite distinct from 
the cases of classical Hermite-Pade approximation known to us. 

Before we proceed to the main body of the paper we need a few preparatory 
lemmas concerning functions f(z) in the complex plane having a Laurent expansion 
around z — oo with finite principal part, 

N 

f(z)= J2 

j=-oo 

For such functions we let II>o and n<o denote the projection operators onto the 
subspaces of nonnegative and negative powers, respectively. For a fixed A we denote 
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by -j^j the Laurent series z^ 1 + Xz~ 2 + X 2 z~ 3 + • • ■. Finally, for Q < I, let lT denote 
the projection onto z~ l . The following two lemmas are simple exercises in the algebra 
of Laurent series. 

Lemma 4.5. Suppose < j and 1 < /. Then 

zJ zJ - A J 



n 



>0 



z — X z — A ' 



z- 



i A-' . z j .V 1 1 



U <°— = — U z-X Z > ■ 

In particular, if p(z) is a polynomial then 

p(z) p{z)-p{\) ,p(z) A'-V(A) 

n ^°^A= z-X ' U z—X = ^^- 

Moreover, 

n0 p(z) _ p(A)-p(0) 
z — A A 
Lemma 4.6. Let k =/= X be given. Then 

, zJ k3+ 1 - 1 - AJ+'- 1 , 

n = z . 

(z — k)(z — A) k — A 

In particular, if p(z) is a polynomial then 

11 (z-«)(z-A) = JTTx Z ■ 

Moreover, 

[p(/s)-p(0) p(A)-p(0) 



n o. P( z ) _ 1 



A 



[z — k)(z — A) k — A 

We will introduce three types of approximation problem, each type specifically 
referring to the pertinent column of a{z) as in Theorem 14.21 The following setup is 
common to all three cases. 

Definition 4.7. Let {Xk,bk, M} be arbitrary numbers satisfying the constraints 

< Ai < A 2 < ••• < A n _i, < 0, M > 0. 

Define (signed) measures 

rl — 1 -. n— 1 

d^(X) = Y,b k S(X-X k ), dv(X) = -—5(X) + Y / Ck6(X-X k ), (4.13) 

fe=i fe=i 
where the c^'s are given by (|3.8p . which amounts to 

- | /(A) diz(A) = M + 1| -M-d M (K)d M (A) (4.14) 
for any function /. Define the two functions 

n-l 



d/i(A) _ ^ 6 



*M = / T=T = E 



z — A ^— ' z — At- 

fc=i 

di/(A) 1/2M Spi c fe 



(4.15) 



z — Afe 

fe=i 



The inverse spectral problem for the discrete cubic string 
and let 



W*{z) = -W(-z) 



MX) 

z + X ' 



Z*(z) = Z(-z) 



di/(A) 
z + A ' 
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(4.16) 



(Note that Z(z) is determined by W(z) and M, and that and Z satisfy (|3.7[) .) 
Finally, let 

K 1 X> 



= / X 3 dfi(X), hj = Iji = 



K + X 



d/i(«)d/i(A). (4.17) 



Definition 4.8. Let A = B = C = V = 1, Ai = I m + tm, 2?i = Po, and 



loo + 2M Iq 



Io,k-l 



B k = 



ho 


hi 


■ ■ ■ h,k- 


ho 


hi 


h.k- 


h-i,o 


Ik-1,1 


■ ■ ■ Ik-l,k- 


ho 


hi ■■■ 


h,k-l 


ho 


hi ■■■ 


h,k-l 


h-i,o 


Ik-1,1 


h-i,k-i 


hi I12 


• ■ ■ hk 




hi I22 


■ ■ ■ hk 





2 < k, 



(4.18) 



1 < k, 



(4.19) 



Ik2 



1 < k, 



(4.20) 



V k = 



ho 


In •• 


■ Ii,fc-i 


ho 


I21 ■• 


• h,k-i 


IkO 


hi ■ ■ 


■ h,k-l 


Po 


ho 


■ • h,k- 




ho 


• • h.k- 


Pk-i 


ho 


■ ■ Ik,k- 



l<k, 



(4.21) 



2 < k. 



(4.22) 



In all these cases, the index k agrees with the size k x k of the determinant. Note also 
that Ak = B k + siyC/i-i for k > 1. 

Theorem 4.9 (Type III approximation problem). Fix 1 < k < n — 1. There are 
unique polynomials Q(z), P(z), P(z) of degree k, k — 1, k — 1, respectively, satisfying 
the following properties: 

(i) Approximation: 



Q 



(ii) Symmetry: 



( w ) ~ ( p ) = 0(z 1} ' as z ^ °°- 

P + PW* + QZ* = 0(z~ (fc+1) ), as z ^ 00. 



(4.23) 
(4.24) 
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The coefficients of Q(z) = 1 + 1j z ^ are given by the k x k linear system 



15 



/ hi 
hi 

V h-i, 



h>2 
hi 

4-1,: 



hk \ 
hk 



Ik-l,k ) 



I qi \ ( ho + w \ 



<12 



'10 



(4.25) 



V ft / 

The polynomials P(z) and P(z) are determined from Q(z) through below. 

Proof. Upon applying n> and Lemma H751 to (|4.23p we obtain 

Q(z)-Q(\) 



P(z)=Tl> (WQ)= J 
P{z)=n> Q (ZQ) = 



z-X 

Q(z)-Q(\) 

z-X 



d/*(A), 
<MA), 



(4.26) 



so that P and P are uniquely determined once we know Q. Inserting this into (|4.24p 
yields 



Tl>o(ZQ) + II>o(l^Q) W* + QZ* = 0(z-( fc+1 >). 



(4.27) 



Applying n> to the left-hand side yields U> ((Z + Z* + WW*)Q), which is zero 
by (|3.7jl . Thus, the nonnegative powers in (|4.27|) vanish automatically, while the 
vanishing of the powers z~ 1 ,...,z~ k is the condition that will determine the k 
unknown coefficients in Q{z). Define the Hankel and Toeplitz operators 

Hf = n <0 o Mf o n> , T f =n> oM/on>o, (4.28) 

where Mf denotes multiplication by the function /. Applying n<o to the left-hand 
side of (|4.27[) yields Hz*Q + Hw*T\yQ, so that Q is determined by the linear system 



{H Z - + H W .T W )Q = 0(z-( fe+1 )). 



(4.29) 



If we restrict the operators to the finite-dimensional subspaces spanned by 



{1, z, . . . , z k } and {z 1 , . 



£ } we obtain k equations 



n'(ff z . + B W .T W )Q = 0, 1 < I < k. 
With the help of Lemma [4751 this can be rewritten as 



n' 



Q(z) 

z + X 



dv(X) 



Q{z)-Q{k) 
(z-k)(z + \) 



d/x(/t)d/x(A) 



0, 



which according to the relation (|4. 14|) between the measures dfj, and du equals 



IT 



Q(z)/z 
2M 



( Q(z)/(z + X) | Q(z)-Q( K ) 
V k + X (z-k)(z + X) 



d/i(/s)d/i(A) = 0. 



(4.30) 
(4.31) 

(4.32) 



Finally, upon executing the projection II' for 1 < I < k using formulas from 
Lemmas 14.51 and 14.01 we obtain 



0(0) 



2M 



goo 

K + x 



d/i(«)d/i(A) = 0, 



(A)'- 1 Q(^) 

K + X 



(4.33) 



d/i(/s)d/i(A) = 0, 



2<l <k. 
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When Q is represented by the column vector (1, g 1; . . . , q^f this takes the form (|4.25|) . 
The determinant of the system (|4.25|) is T>k of Definition 14.81 which is positive by 
Lemma iB.ll in Appendix B. Thus the solution is unique. □ 

Definition 4.10. Denote the results of the fcth approximation problem of type III by 
Qzk , Pzk , Psk ■ 

Corollary 4.11. The highest order coefficient of Qd,k{z) is [Qsk] = (— 1) •Afc/'Cfc- 
Theorem 4.12 (Type II approximation problem). Fix 1 < k < n — 1. There are 
unigue polynomials Q(z), P{z), P(z) of degree k, k — 1, k — 1, respectively, satisfying 
the following properties: 

(i) Approximation: 



Q 



z 



p 
p 



(ii) Symmetry: 

P + PW* + QZ* = 0(z" 
(Hi) Normalization: Q(0) = 0, P(0) = 1. 



0(1) 



as z 



as z 



(4.34) 



(4.35) 



The coefficients of Q(z) = Xw=i Ii 7 -' 1 are given by the k x k linear system 



( ho 

ho 



hi 
hi 



h,k—l 
h.k-1 



\ IkO hi 



Ik.k-1 / 



( 91 \ 

V * J 



Pi 



(4.36) 



The polynomials P{z) and P(z) are determined from Q(z) through l{4-37\ ) below. 

Proof. The proof proceeds in a similar way to the proof of Theorem 14.91 We therefore 
only indicate the steps which are new to the present case. First, from the normalization 
condition we get 

p(z) = 1 + n >0 (WQ) = n>o(iFQ) + 1 - u°(wq), 

Q(z)-Q(X) , (4-37) 



P(z) =n> (ZQ) 

Inserting this into (|4.35|) yields 
n> (ZQ) + n> {WQ) W* + QZ* 



A 



-dv(X). 



(l-U a (WQ))W* =0(2 



-(k+ih 



(4.38) 



Applying n> to the left-hand side yields U> ((Z + Z* + WW*)Q), which is zero 
by (|3.7[) . Thus, the nonnegative powers in (|4.38[) vanish automatically, while the 
vanishing of the powers z~ 1 ,...,z~ k is the condition that will determine the k 
unknown coefficients in Q(z). With the help of the Hankel and Toeplitz operators 
defined in the previous proof we find that Q is determined by the linear system 

{Hz- + Hw*T w )Q + (1 - Tl°(WQ)) W* = 0(z-( fc+1 >), (4.39) 

where we have used that W* = 0(1/ z) according to its definition. If we restrict 
the operators to the finite-dimensional subspaces spanned by {1, z, . . . , z k } and 



: } we obtain the k equations 



If ( (H z * + H W ,T W )Q + (1 - iPiWQ)) W 



0, 



KKk. 



(4.40) 
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With the help of Lemma EQ>1 and the relation (|4.14[) between the measures dfi and dv 
we obtain 



// 



A' Q(k) 
k(k + A) 



dn(K)dfi(X) = / (A^M^A), 



1< I < k. 



(4.41) 



When Q is represented by the column vector (gi, . . . , q^f this gives the system (|4.36j) . 
which has the same determinant T>k > as (|4.25j) . Hence, the solution is again 
unique. □ 

Definition 4.13. Denote the results of the fcth approximation problem of type II by 
Qzk+i , Psk+i , P3k+1 ■ 

Corollary 4.14. The highest order coefficient of Q3k+i(z) is [Qsk+i] = 
{-if-^V'jVu. 

Theorem 4.15 (Type I approximation problem). Let < k < n— 1. Then there are 
unigue polynomials Q(z), P{z), P(z) of degree k + 1, k, k, respectively, satisfying the 
following properties: 

(i) Approximation: 

0(1) 



Q 



z 
w 



p 
p 



as z 



(ii) Symmetry: 

P + PW* +QZ* = 0(z- (fc+1) ), asz 
(Hi) Normalization: Q(0) = 0, P(0) = 0, P(0) = 1. 



(4.42) 
(4.43) 



ofQ(z) = E' 


=1 % 7j 


are given by the {k - 


-1) 


x (k + 1) linear system 


( J oo + m 


hi ■ 


•• hk \ 




( ql ) 






( 1 ^ 




ho 


hi ■ 


■ ■ hk 















I20 


hi ■ 


■ ■ hk 













(4.44) 


\ ho 


hi ■ 


■ ■ hk ) 




\ Qk+i J 






V j 





The polynomials P(z) and P(z) are determined from Q(z) through fr4-45\ ) below. 
Proof. From the normalization and approximation conditions we get 

p(z) = n >0 (wQ) =n> (^g) - n°(WQ), 

Q(z)-Q(X) , „ s (4.45) 



z-X 



-dv(X). 



P(z)=n> Q (ZQ) = 
Inserting this into (|4.43|) yields 

n> a (zQ) + n> (i^g) w* + qz* - ii (wq) w* = o(z 



-(fe+i; 



(4.46) 



Applying n> to the left-hand side yields H> ((Z + Z* + WW*)Q), which is zero 
by (|3.7[) . Thus, the nonnegative powers in (|4.46p vanish automatically, while the 
vanishing of the powers . . . , z~ k will give k conditions on the k + 1 unknown 
coefficients in Q(z), which together with the normalization condition P(0) = 1 is the 
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required number of conditions. Symbolically, with the help of the Hankel and Toeplitz 
operators defined earlier, we find that Q is determined by the linear system 

(H z * + H W *T W )Q - U°(WQ) W* = 0( Z -^), 

n°(ZQ) = 1. 

If we restrict the operators to the finite-dimensional subspaces spanned by 
{z, . . . , z k+1 } and {z -1 , . . . , z~ k }, then the k equations alluded to earlier are 

\(H Z * + H W ,T W )Q - U°(WQ) W*) = 0, 1 < I < k. (4.48) 



Likewise, Lemma 14.61 and the relation (|4.14p between the measures d/i and dv imply 
that the complete set of equations determining Q reads 

Q ^ M x) = i, 

(4-49) 



A 

X 1 Q(k 



-d/i(«)d/i(A) = 0, 1 < I < k. 



k(k + A) 

When Q is represented by the column vector (qi, . . . , qk+i)* this gives (I4.44|) . The 
determinant of the matrix of coefficients is Ak+i = Sfe+i + 2lfCk in the notation of 
Definition 14.81 Following the methods used in 22, Lemma 4.10] one can establish 
Heine-type formulas for both determinants in this sum, which show that both are 
strictly positive provided k < n — 1; see Lemmas IB. 2 1 and IB. 31 in Appendix B. From 
these formulas one also concludes that when fc = n — 1 the first determinant is zero, 
since the measure d/i(A) is supported onn-1 points, while the second remains strictly 
positive. This proves that the system (|4.44[) has a unique solution for the claimed range 
of k. □ 

Definition 4.16. Denote the results of the fcth approximation problem of type I by 
Q3/C+2 , -P3/C+2 , -P3/C+2 ■ 

Corollary 4.17. The highest order coefficient of Qzk+2{ z ) is [Q^k+2\ = 
(~l) k+1 C k /A k . 

The last approximation problem of type I has a particular significance. 
Corollary 4.18. For k = n — I, the approximation problem of type I gives 

z 



V3n~l 

Moreover, the approximation is exact: 



(z) = -2Mz\{ (l- -). (4.50) 



z{z) = W{z) = (4.51) 

Q3n-l(Z) Q3n-l(Z) 



Proof. We write the conditions (|4.49p that Qs n -i has to satisfy: 

Q3n-l(A) 



A 

X Q3n-l(«) 



k(k + A) 



di/(A) = 1, 

d//(«)d/i(A) = 0, l<«<n-l. 



(4.52) 
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We claim that Q3„_i has to vanish on the support of the measure d/x. We recall 
that Qsn-i is a polynomial of degree n with one root at A = by the definition 
of the approximation problem of type I. Thus the claim is that the remaining roots 
are precisely the Xj occuring in (|4.13p . To prove that, we write explicitly the second 
condition, 



n-l 

^ V./ A, "'""' ) V Aj 



Xl - M x)) ( 9^Mh i =0 , i < i < n-i. 



It is straightforward to establish a multiple integral formula of Heine type for the 
matrix of coefficients, from which one concludes that the matrix is nonsingular; see 
(IB . 13|) in Appendix B. This implies that Qzn-i{Xj) = for 1 < j < n — 1, which 
proves the claim. Moreover, the definition of the measure di' and the hrst condition 
in (|4.52|) imply that the coefficient of the linear term of Q3 n -i is — 2M, which proves 
the product formula for Qa n -i. To prove that the k = n — 1 approximation of type I 
is exact we show that 

n <0 (Wxf % n - 1 ))=o. 



W J V P3n-1 

This follows from Lemma 14.51 since for example 



n <0 (Q 3n -iW)( z ) = 1 9^M dKX ) = J Q3 ;:^ A) d M (A) = o 

because Q3 n ~i vanishes on the support of dfi. Likewise, 

n°(Q 3n - lW ) =no/ ^^d M (A) = / ^iHdM(A) = 0. 

The rest follows from the definition of the approximation problem of type I. □ 

We can now present the solution of the uniqueness part of the inverse problem. 
Theorem 4.19. Let 

< Ai < A 2 < ••• < A„_i, &!,...,&„_! <0, M> 0(4.53) 

be the spectral data (as described in Section [3|) of a cubic string {m^, Xfc}^ =1 with 
mk > and x\ < ■ ■ ■ < x n . Then the string data can be recovered uniquely from the 
spectral data, up to a translation along the x axis, from the determinantal formulas 

CkVk ln- k = (4.54) 



2A k+1 A k ' V' k 
where lj = Xj+i — Xj. (See Definitions \4- 7| and \4--8\ for notation.) 

Proof. Equations (|4.54[) follow from Corollaries 14.111 14. 441 14.171 and equations (|4.1ip 
and (1442)1 . □ 

To show the existence part (that any numbers satisfying the constraints ()4.53l) 
really are spectral data of some cubic string), we need the following theorem, whose 
proof will be published elsewhere. Note that the orthogonal polynomials occurring in 
the description of continued fractions of Stieltjes type satisfy three-term recurrence 
relations. It is therefore gratifying to know that the polynomials occurring in the 
approximation problems for the cubic string satisfy four-term recurrence relations. 
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Theorem 4.20. Consider the same setup as in Definition ^. 7| (numbers {X k ,b k , M} 
given), and let the polynomials Qj(z), Pj(z), Pj(z) be defined (for j > 2) by the 
unique solutions of the approximation problems of type I, II and III, as in Definitions 
\4-10\ and \l.lb] Define lj and rrij from the numbers {X k ,b k , M} via ft4.54\ ), and 
consider the recurrence relations 

X3k — 2^n-k ^3fc-l + ln-k ^3/c-2 + ^3fc-3, 1 < < 71 — 1, 

X 3k+X =l n - k X 3k -i+X 3k -2, l<k<n-l, (4.55) 

X 3k+2 = -2zm n - k X 3k +X 3k - 1 , 0<fc<n-l. 

Then with initial conditions (X-i,Xo,X\) = (1,0,0) the solution (for j > 2) is 
Xj = Pj , when (X—i,Xa,Xi) = (0,1,0) the solution is Xj — Qj, and when 
(X_i,X ,Xi) = (0,0, 1) the solution is X.j = Pj. 

Theorem 4.21. Given any numbers {X k ,b k , M} satisfying the constraints ^4.53\ ), the 
formulas fr4-54\ ) define a cubic string with these numbers as spectral data. 

Proof. Note first that all rrij and lj denned by (|4.54[) are strictly positive. This follows 
since V' k < according to (|B.5|) and our definition of d/z(A) as a negative measure, 
while the remaining determinants are strictly positive. This means that we lrave a valid 
cubic string, but we do not know yet whether its spectral data {X k ,b k , M} coincide 
with the numbers {Xk,bk, M} that we started with. However, plugging {X k ,b k , M} 
into (|4.54p recovers the same string data {wij,lj} by uniqueness (Theorem 14. 1 9() . 
Thus, applying Theorem 14.201 to the set of numbers {X k ,b k , M} produces the same 
polynomials Qj, Pj, Pj as one would get with {X k , b k , M}, since the string data are the 
same in both cases. And according to Corollarv l4.18l both M and W(z) are uniquely 
determined from these polynomials. Thus = X k , b k = b k , and M = M, which is 
what we needed to show. □ 
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Appendix A. A collection of facts about totally nonnegative and 
oscillatory matrices 

We proved that the spectrum of the Neumann-like discrete cubic string (|l.ip is positive 
and simple by rewriting it as a matrix eigenvalue problem with an oscillatory matrix. 
Below, the reader will find some relevant definitions and theorems from the vast area 
of totally nonnegative matrices; for details see the book by Gantmacher and Krein 
|26j and the excellent survey by Fomin and Zelevinsky |27j . 

Definition A.l. A square matrix is totally nonnegative (TN) if all its minors are 
nonnegative. It is totally positive (TP) if all its minors are positive. It is oscillatory 
if it is TN and some power of it is TP. 
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Theorem A. 2. All eigenvalues of a TP matrix are positive and of algebraic 
multiplicity one, and likewise for oscillatory matrices. All eigenvalues of a TN matrix 
are nonnegative, but in general of arbitrary multiplicity. 

Theorem A. 3. The product of an oscillatory matrix and a nonsingular TN matrix 
is oscillatory. 

Theorem A. 4 (Star operation). If A = (djj) is oscillatory, then (A*)^ 1 is oscillatory, 
where A* = a^) . 

Theorem A. 5. A symmetric tridiagonal matrix is oscillatory if and only if it is 
positive definite and has positive off- diagonal entries. 

Definition A. 6. A planar network (T,ui) of order n is an acyclic, planar directed 
graph r with arrows going from left to right, with n sources (vertices with outgoing 
arrows only) on the left side, and with n sinks (vertices with incoming arrows only) 
on the right side. The sources and sinks are numbered 1 to n from top to bottom. 
All other vertices have at least one arrow coming in and at least one arrow going out. 
Each edge e of the graph T is assigned a scalar weight uo(e). 

Definition A. 7. Given a planar network T of order n, the weight of a directed path 
in r is the product of all the weights of the edges of that path. The weighted path 
matrix A{T, lu) is the n x n matrix whose entry is the sum of the weights of the 
possible paths from source i to sink j. 

The following famous theorem was discovered independently by Lindstrom, 
Karlin-McGregor, and Gessel-Viennot. 

Theorem A. 8. Let I and J be subsets of {1, . . . ,n} with the same cardinality. The 



of the weighted path matrix A(T, uj) of a planar network is equal to the sum of the 
weights of all possible families of vertex- disjoint paths connecting the sources labelled 
by I to the sinks labelled by J . (The weight of a family of paths is defined as the 
product of the weights of the individual paths.) 

Corollary A. 9. // all weights u) are nonnegative, then the weighted path matrix is 



minor 



Ai,j(A) = det(ay) ie /^- e j 



TN. 



Appendix B. Heine-like determinant formulas 



The following was proved in [55], where it appeared as Lemma 4.10. 
Lemma B.l. Suppose /i is a measure on R such that the integrals 




(B.l) 



are finite. Let 




(B.2) 
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where 

A(x) = A(xi, ...,x k ) = ^Q(x, - Xj), 

(B.3) 

T(x) =T(xi,...,Xk) = Y\(xi+Xj). 

i<j 

(When k = or 1, we let A(x) = T(x) = 1. Likewise, uq = vq = 1 .) Then for k > 1 
£/ie following k x k determinant formulas hold: 
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The following lemmas state some additional determinant formulas needed in this 
paper. See Definition 14.81 for notation. 



Lemma B.2. For k > 1, 



where 



B, 



A/ = A(x n , 




^A^iyrz^ dfi 2k (x), (b.7) 



; • • • 3 ^ik )j 



where the sum runs over all ( 2 ^ fc ) k-element subsets I — {ii, . . . ,i k } of {1, . 
and where I c denotes the complement of I. Moreover, letting 

A(x) 2 dfi k (x) 



1 

h — TT 



k\ J R k T(x) X±X2---Xk 



for k > 1, 



we have 



B k = 



2 k 



t k +i 



Uk-l 
Uk 



(B.8) 

..,2k}, 

(B.9) 
(B.10) 



The proof of (|B.7|) is omitted since it is very similar to the proofs in . Equation 
(|B.10p follows by rewriting the integrand in (|B.7|) using the symmetric function identity 
(2.76) from [22] . By simply replacing the measure dfi(x) by id/i(i) we obtain the 
corresponding formulas for C k : 
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Lemma B.3. For k > 1, 
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Moreover, 



C k = 
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(B.11) 



(B.12) 



It follows from (|B.7|) and (|B. 1 1|1 that if dfi is a discrete measure supported on 
n — 1 positive points, like in the main text, then Bk and Ck are strictly positive for 
1 < k < n and zero for k > n (because then each A/ vanishes). 

Finally we state a formula for the determinant appearing in the proof of 
Corollary 14.181 We omit this proof as well; it is also similar to the previous ones 
(and actually a little easier, since no double integrals are involved). 

Lemma B.4. Let the measure dfi be the same as in the main text and define the 
(n — 1) x (n — 1) determinant £ = det[£y ] where 



x + A 



Then 



A(A n -l, A n -2, • • ■ , Ai) 
(n-1)! 



A(x) 2 



R"" 1 Hi,j=i(vi + Aj) 



d^"- 1 ^). (B.13) 



In particular, £ ^ 0. 
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